Cavendish-HEP-2004/1 1 

DAMTP-2004-26 

DESY 04-085 

The Gluon Green's Function in N = 4 
Supersymmetric Yang— Mills Theory 



'^ ; Jeppe R. Andersen ^'^ , Agustm Sabio Vera '^ ^ 



O 






X 



Cavendish Laboratory, University of Cambridge, Madingley Road, CBS OHE, 

Cambridge, UK 



.<^ ; ^DAMTP, Centre for Mathematical Sciences, Wilberforce Road, CBS OWA 

Cambridge, UK 



"^II. Institut fiir Theoretische Physik, Universitdt Hamburg, Lumper Chaussee 149, 
CN ■ 22761 Hamburg, Germany 

a^ 
o 
o 

^ ■ Abstract 

o 

(<— ^ I The high energy Hmit of scattering amphtudes in the N = 4 supersymmetric Yang- 

Mills theory is studied by solving the corresponding BFKL equation in the next- 
to-leading approximation. The gluon Green's function is analysed using a newly 
CL|| proposed method suitable for investigating the contribution from higher conformal 

^ ' spins. From this new approach complete agreement is obtained with the results of 

Kotikov and Lipatov on conformal spins and angular dependence. 



1 Introduction 



In recent years there has been a large interest in the study of the next-to- 
leading (NLL) corrections to the Balitsky-Fadin-Kuraev-Lipatov [1] (BFKL) 
equation in QCD [2-5]. The calculation of the NLL corrections was extended 
to supersymmetric gauge theories in Ref. [6]. In that work it was shown how 
in a N = 4 SYM theoretical playground the kernel of the BFKL equation 
is simplified. In particular, the analyticity of the eigenvalues in terms of the 
conformal spins allowed the study of the connection between the DGLAP and 
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BFKL equations in this model, this being possible mainly due to the lack of 
coupling constant renormalisation. In Ref. [7] the anomalous dimension matrix 
of the Wilson twist-2 operators in this maximally supersymmetric theory in 
four dimensions was calculated at two loops. Very recently these results were 
extended to three loops in Ref. [8] using the calculation for the non-singlet 
case first obtained in Ref. [9], and in agreement, in the supersymmetric limit, 
with the results for the singlet case derived in Ref. [10]. 

In the present study, the method for solving the BFKL equation at NLL ac- 
curacy developed in Ref. [14, 15] is applied to N = 4 SYM. This will serve two 
purposes. Firstly, since the method solves the BFKL equation with full angu- 
lar information, it can be used to test the results in the literature [6] for the 
dependence of the eigenvalues of the kernel on conformal spins. Secondly, due 
to the conformal invariance of the theory also at higher orders, the analytic so- 
lution to the NLL BFKL equation is known for N = 4 SYM, and therefore it is 
possible to directly test the results obtained using the completely independent 
method described in this work. The results in the literature for the solution of 
the N = 4 SYM BFKL equation at NLL in terms of conformal eigenf unctions 
are reviewed in the next section. In Sec. 3 the NLL BFKL equation is solved 
directly in momentum space using the iterative method of Ref. [14,15]. Sec. 4 
is devoted to a study of the numerical structure of the trajectory and the real 
emission kernel arising in this approach. In Sec. 5 the results obtained in the 
two approaches are compared and the conclusions are presented. 



2 The BFKL equation and its solution in Mellin space 



In the BFKL formalism the high energy limit of a scattering process factorizes 
as 



^(^)= /^ /^ '^(^'^) ^'ih)f(ka,h,Y^\n^ 

where sq = |A:a||^fe| is the Regge scale. The energy dependence is determined 
by the universal process-independent gluon Green's function /. The impact 
factors, $, $', depend on the process under study. In the Regge-limit of large 
centre of mass energy and fixed momentum transfer the effective degrees of 
freedom are the transverse momenta ka^b of the exchanged gluons. The dynam- 
ics of these processes is then two-dimensional evolving with a time variable Y. 
The evolution of the gluon Green's function with Y is governed by the BFKL 
equation. This equation is traditionally written in terms of a Mellin transform 
of Y, i.e. 
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With such transformation the equation reads 
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The inhomogeneous term of this integral equation corresponds to a single 
gluon exchange and the kernel (in dimensional regularization, Z^ = 4 + 2e) 
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contains the gluon Regge trajectory [6], 
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with the coupling a = ^^r, which is not running in N = 4 SYM. The kernel 
also includes two contributions from the real emissions: 
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In N = 4 SYM the BFKL kernel respects conformal symmetry even at NLL 
accuracy, and the eigenfunctions do not change compared to the leading loga- 
rithmic (LL) ones. The solution to this equation can therefore be found using 
the expansion on the known eigenfunctions 
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with 6 the angle defined by the ka and kb transverse momenta. The Fourier 
transform in angles is characterized by the so-called conformal spins n. Given 
that the coupling is fixed it is possible to fully diagonalize the BFKL kernel 
by simply acting on the LL eigenfunctions, i.e., in the MS scheme we have [6] 



d'qlC{k,q)(^) e^"^ = cu„(a,7)=e'''x(l^|,7) + r7 + fi(|n|,7).(10) 



In this expression we have defined 
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We have checked that switching to the gluon-bremsstrahlung (GB) scheme 



aGB = a+(^-^ja\ ^''^ = aGB, (16) 

or to the dimensional reduction (DRED) scheme, which respects SUSY [6], 



«DRED — fl+— a , C — «DRED i^OORED , (17) 

the results of this work do not change qualitatively. 

The analyticity of these expressions, obtained by Kotikov and Lipatov in [6], 
for the dependence of the eigenvalues of the NLL BFKL kernel on the confor- 
mal spins is very important. It allows the analytic continuation to negative \n\ 
and to find the connection between the BFKL equation and DGLAP in the 
N = 4 supersymmetric gauge theory. One of the objectives of the present work 
will be to confirm that these results for the conformal spins are correct using 
a completely orthogonal method of solution of the BFKL equation, which is 
developed in the next section. 



3 The solution directly in transverse— momentum space 



The starting point of the new method of solution is Eq. (3) with the infrared di- 
vergences regularized in dimensional regularization. In the real emission kernel 
there are contributions which will lead to e poles after phase space integration, 
}Cl!'\ and others which will be finite, JCr, i.e. 
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In order to explicitly show the cancellation of the infrared e divergencies we 
introduce a phase space slicing parameter. A, in the integral over real emission, 
and make use of the approximation 



/. [k + ka, h) -f^[k + ka, h) e [e - \^) + /^ (A;,, h) e (A^ - e) . (19) 
This is valid when A is small compared to ka- Hence the equation now reads 
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C0f^{ka,h)=5^'+''^{ka-h) (20) 
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We have performed the integration over phase space for those emissions below 
the infrared cut-off A, obtaining the result 



ld^+HlCi'^{q,q + k)d{x''-P 



a 1 /\^' 



^{^"KS 



--2C(2) + e(^-- + 14C(3) 



(21) 



Now it is possible to show how the poles cancel by calculating the gluon Regge 
trajectory in our regularisation scheme, i.e. 
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where ^ and rj coincide with those defined in Eq. (11). The corresponding real 
emission part will be of the form 
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Finally the N=4 SYM NLL BFKL equation can be expressed as 



(cj - LOO (Pa, A')) /. (ka, 4) = 6^'^ (ka " h) 

dPk (^^ e""' (P - A^) + iCr (ka, ka + k)) fu (K + K h) . (24) 

Note that the treatment of the kernel is with its full angular dependence, i.e. 
without angular averaging over the angle between ka and A;^. It is therefore 
possible to extract the contribution to the solution from all conformal spins. 

Eq. (24) will be solved by iteration, generalising the procedure of Ref. [11-13]. 
In order to do so it is useful to introduce the notation 
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(25) 



The uj dependence can go to the denominator of the right hand side of the 
equation and iterate, generating in this way multiple poles in the complex u 
space, i.e. 
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It is now possible to invert the Mellin transform as in Eq. (2) and go back to 
energy space, with the final compact expression for the gluon Green's function 
being 
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where the notation t/o = ^ has been used. We note that the solution has been 
expressed as the phase space integral of a product of effective emission vertices 
connected with no-emission probabilities. 

Before proceeding further in the numerical study of the solution we present 
an analysis of the trajectory and real emission kernel in the next section. 
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Fig. 1. Comparison of the LL and NLL gluon Regge trajectory calculated in different 
schemes. 

4 Analysis of the trajectory and emission kernel 



In the following we indicate which expressions have been used in the imple- 
mentation of the solution to the BFKL equation. Firstly, in Fig. 1 it is shown 
the behaviour of the gluon Regge trajectory as in Eq. (22), ujo{(f,X), as a 
function of A and q. It is interesting to note that the LL trajectory always lies 
below the NLL one, for all regular isat ion schemes. This is the opposite effect 
to that found in the QCD case, see Ref. [15]. The correction to the trajectory 
is smallest in the GB scheme, with the DRED and MS being very similar to 
each other. For the A dependence we can see that, at a fixed value of g = 20 
GeV, the negative value of the trajectory decreases for lower values of A. Due 
to the logarithm, the trajectory also decreases when, for a fixed A = 1 GeV, 
q increases. For these plots the coupling was chosen to be a = 0.2 for all 
schemes. The plots at the right hand side of Fig. 1 show the ratio of the NLL 
uJo{(f, A) to its LL value. 

Secondly, the part of the real emission kernel in Eq. (8), can be written as 
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with 9 being the angle between the two-dimensional vectors g and if. For the 
function JF (g, q', 9) we use the expression 
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As in the QCD case, see Ref. [15], this kernel has integrable singularities at 
q = (f ^q = and g* = 0. This structure is revealed when the kernel is plotted as 
in Fig. 2 where /C^ (g, g' = 20 GeV, ^) is shown for a = 0.2 in the MS scheme. 
Note how this kernel is positive in a large region of phase space, contrary to 
the QCD case [15]. 



5 Study of the gluon Green's function 



The eigenvalues of the NLL BFKL kernel as in Eq. (10), cj„(a, 7), along the line 
7 = \ + ii' for a value of the coupling of a = 0.2 are plotted in Fig. 3. This line 
parameterises the contour of the 7 integration in Eq. (9). The LL eigenvalues 
are compared to those obtained at NLL for several values of the conformal 
spin. At high energies, for n > the relevant region is the one close to z/ = 0. 
The zero conformal spin evolution is governed by the two maxima and is the 
dominant contribution among all conformal spins. We will come back to this 
figure when the evolution with energy of the different contributions to the 
gluon Green's function is studied below. 

Very importantly, for a Regge-like choice of energy scale, the NLL kernel 
in the N = 4 SYM theory does not develop an imaginary part along the 
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Fig. 2. Structure of the kernel K^ {q, q', 9) for q' = 20 GeV as a function of q and 
the angle between g and q' for the coupling a = 0.2. 




Fig. 3. The eigenvalues of the BFKL kernel as a function of 7 = g + i^ at LL and 
NLL for different values of the conformal spin n = 0,1,2,3,4 and a coupling of 
a = 0.2. 
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Fig. 4. The gluon Green's function for zero conformal spin for fixed values of the 
external scales, a coupling of a = 0.2 and large values of the energy scale Y. 

7 = i + zz^ contour. This implies that the asymptotic behaviour at large Y is 
well controlled and that the gluon Green's function is monotonically growing 
with energy. To illustrate this point, we plot in Fig. 4 the angular averaged 
gluon Green's function up to very large Y both at LL and NLL accuracy. 

With the intention to confirm the calculation of the conformal spins in N = 4 
SYM of Ref. [6] the two methods of solution, that in Mellin space described in 
Section 2, and the one proposed in this work directly in transverse momentum 
space as explained in Section 3, will be used to calculate the behaviour of the 
Green's function for different conformal spins. To proceed, Eq. (9) can be 
written as 






(31) 



and the coefficients in the expansion can be calculated following the solution 
in Mellin space, i.e. 
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or they can be obtained using the solution in k space, for this it is needed to 
project on angles, i.e. 
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Fig. 5. Contributions to the gluon Green's function from different conformal spins 
as a function of energy. 



/n(|4|,|4|,Y) = j — f (ka,kb,Y) cos{n9). 



(33) 



Therefore the projections on conformal spins can be compared to each other 
using two completely independent approaches. The results obtained from both 
solutions are shown in Fig. 5. Both methods exactly coincide in their predic- 
tions. This confirms the validity of the results calculated in Ref. [6] and it is 
a very serious test of the method first proposed in Ref. [14, 15] to study the 
gluon Green's function at NLL. As expected from Fig. 3, the dominant confor- 
mal spin is n = 0, whose corresponding eigenvalue is the only positive one at 
z/ = in Fig. 3. The eigenvalue at z/ = for the n = 1 conformal spin is zero, 
and thus it is expected to give a constant contribution at large energies. This 
behaviour is indeed observed in Fig. 5. For the rest of conformal spins their 
contributions to the gluon Green's function decrease as the available energy 
is larger, their eigenvalue being negative at the vecinity of z^ = 0. 

A natural question is the convergence in n of the angular expansion of Eq. (31). 
As the method of solution proposed in this paper allows for a full determina- 
tion of the angular dependence it is possible to answer this point in a simple 
manner. This issue is addressed in Fig. 6 where the gluon Green's function is 
plotted as function of the angle between the two transverse momenta. Here it 
can be seen that the conformal expansion reaches good convergence for con- 
formal spins above n = 8 for the chosen values of Y. This plot confirms again 
that both approaches produce exactly the same results. The graph is produced 
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Fig. 6. The dependence of the gluon Green's function on the angle between the 
transverse momenta ka and kb- 

for two different energies showing a stronger angular correlation for lower en- 
ergies (the curve is flatter in 9 for Y = 5), a consequence of the increasing 
dominance of the zero conformal spin at larger energies. 



6 Conclusions 



The solution to the NLL BFKL equation for the N = 4 supersymmetric Yang- 
Mills field theory has been studied in detail. In particular, a newly proposed 
method of solution has been used which allows for a detailed study of the 
dependence of the gluon Green's function on the conformal spins. The results 
of this work confirm that the calculations of Ref. [6] are correct and, simul- 
taneously, that the proposed method of solution of Ref. [14, 15] for the NLL 
BFKL equation provides the true answer and accurate description of angular 
dependences in the multigluon ladder. The growth with energy of the gluon 
Green's function for the Regge-like choice of scale has also been demonstrated. 

Acknowledgements We thank Tolya Kotikov and Lev Lipatov for very useful 
discussions, the CERN Theory Division for hospitality, and the IPPP, Univer- 
sity of Durham, for use of computer resources. J.R.A. wishes to thank the II. 
Institut fiir Theoretische Physik, University of Hamburg, and A.S.V. thanks 
the Cavendish Laboratory at the University of Cambridge for hospitality. 



13 



References 



[1] L. N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338 [Yad. Fiz. 23 (1976) 642], 
E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 45 (1977) 199 
[Zh. Eksp. Teor. Fiz. 72 (1977) 377], 

I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822 [Yad. Fiz. 
28 (1978) 1597]. 

[2] L.N. Lipatov, JETP 63, 904 (1986), 

G. Camici and M. Ciafaloni, Phys. Lett. B 395, 118 (1997), 
R. S. Thome, Phys. Lett. B 474 (2000) 372, Phys. Rev. D 64 (2001) 074005, 
J. R. Forshaw, D. A. Ross and A. Sabio Vera, Phys. Lett. B 498 (2001) 149, 
M. Ciafaloni, M. Taiuti and A. H. Mueller, Nucl. Phys. B 616 (2001) 349, 
M. Ciafaloni, D. Colferai, G. P. Salam and A. M. Stasto, Phys. Lett. B 541 

(2002) 314, Phys. Rev. D 66 (2002) 054014. 

[3] V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429 (1998) 127. 

[4] M. Ciafaloni and G. Camici, Phys. Lett. B 430 (1998) 349. 

[5] D.A. Ross, Phys. Lett. B431 (1998) 161, 

Yu.V. Kovchegov and A.H. Mueller, Phys. Lett. B439 (1998) 423, 

J. Bliimlein, V. Ravindran, W.L. van Neerven and A. Vogt, preprint DESY-98- 

036, hep-ph/ 9806368, 

E.M. Levin, preprint TAUP 2501-98, hep-ph/9806228, 

N. Armesto, J. Bartels, M.A. Braun, Phys. Lett. B442 (1998) 459, 

G.P. Salam, JHEP 8907 (1998) 19, 

M. Ciafaloni and D. Colferai, Phys. Lett. B452 (1999) 372, 

M. Ciafaloni, D. Colferai and G.P. Salam, Phys. Rev. D60 (1999) 114036, 

R.S. Thorne, Phys. Rev. D60 (1999) 054031, 

S. J. Brodsky, V. S. Fadin, V. T. Kim, L. N. Lipatov and G. B. Pivovarov, 

JETP Lett. 70 (1999) 155, 

C. R. Schmidt, Phys. Rev. D 60 (1999) 074003, 

J. R. Forshaw, D. A. Ross and A. Sabio Vera, Phys. Lett. B 455 (1999) 273, 

G. AltarelH, R. D. Bah and S. Forte, Nucl. Phys. B 575 (2000) 313, Nucl. Phys. 

B 621 (2002) 359, Nucl. Phys. B 674 (2003) 459, 

M. Ciafaloni, D. Colferai, G. P. Salam and A. M. Stasto, Phys. Lett. B 576 

(2003) 143, Phys. Rev. D 68 (2003) 114003, Phys. Lett. B 587 (2004) 87. 

[6] A. V. Kotikov and L. N. Lipatov, Nucl. Phys. B 582 (2000) 19, Nucl. Phys. B 
661 (2003) 19 [Erratum-ibid. B 685 (2004) 405]. 

[7] A. V. Kotikov, L. N. Lipatov and V. N. Velizhanin, Phys. Lett. B 557 (2003) 
114. 

[8] A. V. Kotikov, L. N. Lipatov, A. I. Onishchenko and V. N. Velizhanin, 
hep-th/0404092. 

[9] S. Moch, J. A. M. Vermaseren and A. Vogt, Nucl. Phys. B 688 (2004) 101. 



14 



[10] A. Vogt, S. Mocli and J. A. M. Vermaseren, |hep-ph/0404111| 



[11] J. Kwiecinski, C. A. M. Lewis and A. D. Martin, Phys. Rev. D 54 (1996) 6664 

[12] C. R. Schmidt, Phys. Rev. Lett. 78 (1997) 4531 

[13] L. H. Orr and W. J. Stirhng, Phys. Rev. D 56 (1997) 5875 

[14] J. R. Andersen and A. Sabio Vera, Phys. Lett. B 567 (2003) 116. 

[15] J. R. Andersen and A. Sabio Vera, Nucl. Phys. B 679 (2004) 345. 



15 



